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We (analytically) calculate the energy spectrum corresponding to various exper- 
imental and numerical turbulence data analyzed by Benzi et al. [Q. We find two 
bottleneck phenomena: While the local scaling exponent Cr(^) of the structure func- 
tion decreases monotonically, the local scaling exponent ( p {p) of the corresponding 
spectrum has a minimum of CpiPmin) ~ 0.45 at p m in ~ (lOr?) -1 and a maximum of 
CpiPmax) ~ 0.77 at Pmax ~ 8L _1 . A physical argument starting from the constant 
energy flux in p-space reveals the general mechanism underlying the energy pileups 
at both ends of the p-space scaling range. In the case studied here, they are induced 
by viscous dissipation and the reduced spectral strength on the scale of the system 
size, respectively. 
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In nonlinear dynamics scaling exponents in r-space and p-space are often identi- 
fied with each other. In the case of fully developed turbulence the scaling exponent 
£r,2 of the velocity structure function D^ 2 '(r) = ((u(x + r) — u(x)) 2 ) oc r^* 2 is 
believed to coincide with the scaling exponent ( Pt2 of the energy spectrum of the 
velocity field (multiplied by p), E{p) oc p^^ 2 ^ 1 . In this letter we demonstrate 
that this identification only holds in the limit of very large Reynolds number. For 
(Taylor)-Reynolds numbers Re\ < 200 typically achieved in full numerical simu- 
lations bottleneck phenomena || lead to considerable differences between 
and Cp,2 • First, for large p near the (inverse) scale of dissipation, the bottleneck 
effect accounts for the puzzling observation that the numerical spectrum is clearly 
flatter both in experiment [f|, || and in numerics 0, |6| than E(p) oc p~ 5 / 3 instead of 
being steeper as expected from the possibility of intermittency corrections. Second, 
for small p near the inverse external length scale L _1 (where L is defined by the 
driving force), a similar bottleneck effect leads to steeper spectra. This is another 
hint that finite size effects as also found and analyzed in |7|, |8], |9], [10[] have to be 



considered. Our observation has far-reaching consequences for both the numerical 
and experimental determination of asymptotic scaling exponents from spectra. 

We first focus on the crossover between viscous subrange (VSR) and inertial 
subrange (ISR) and start from Benzi et al.'s |l| measured longitudinal |[l l] velocity 
structure function, assuming that the system size L — > oo M, i.e., no large scale finite 
size effects [[F[ [L2| are considered. Benzi et al. Jl| analyzed various numerical and 
experimental data by means of the extended self similarity method [[13], |l| and found 
that for r > rj (77 being the Kolmogorov scale) the m th longitudinal velocity structure 
function D^fr) obeys D^ir) = C m (r/(r/^)) < ' r ' m , with a universal function f(r/rj) 
for all moments m, for all Re\, and for all kinds of isotropic flow. We restrict 
ourselves to the second order structure functions and drop the index 2 in what 
follows. 

The structure functions -D^(r) and D(r) are connected by |TT| D(r) = 3Di(r) + 



DL(r)d\nDL(r)/d\nr. Both functions can be fitted by parametrizations of the 
Batchelor type JT3], O. Originally given by Batchelor as a parametrization, this 
formula (1) recently got theoretical support by Sirovich, Smith, and Yakhot 0, 
who moreover find agreement between the Batchelor energy spectrum and numeri- 



cal spectra [16, for 30 orders of magnitude. The high quality of the Batchelor fit 
has also been established by older experiments, for an overview see [jTI|, 13 1. Here, 
the Batchelor fit for D(r), 



D(r) r 2 /(3r] 2 



V 2 A /lX 3/ 2/ ,2\l-C/2 : 



3b J \rj 



is slightly superior to that of Di(r), see fig. 1. Here, ( denotes the asymptotic 
value of ( r for r >i), and v v and b are the Kolmogorov velocity and Kolmogorov 
constant |TT[| , respectively. The experimental value f(r — rj) = 8.577 • 10~" 



f(l) determines b = ll/(45(/(l)) 2/3 ) = 5.834, slightly smaller than b w 6.0 - 8.4 
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found in older experiments, which also show excellent agreement with ([[]) [11], [14 . 



We determine D(r) from a spline fit to the data and compare the result with the 
Batchelor parametrization ([I]) in fig.l. There are no visible deviations. 
The /ocaHogarithmic slope [0] of eq. (|TJ), 



rflnP(r) (2-Cy 
UJ ~ d\nr ~ r 2 + r 2 



is monotonically decreasing for increasing r. Here, r d = (36) 3/4 r/ (for ( = 2/3) 
determines the r-space crossover, defined by equating the limits for large and small 
rofeq. ©, (r d / V ) 2 /3 = b(r d / V ) 2 / 3 . 

Next, we calculate the spectrum E(p) which is, when neglecting boundary terms, 
given by |ll], [12 

1 r°° 

E(p) = / prsin(pr)D(r)dr. (3) 

2% Jo 

In view of our results in fig.l we feel justified to consider eq. ([I]) as an exact descrip- 
tion of the experimental structure function of ref.fllj. Inserting eq. (|l|) into eq. (0) 
we obtain 

, ■ d 3 r°° exp(ipr d x) 

[P) ~ 12tt V 2 dp 3 J ~oc (i + x 2 y-v 2 

r 3o-l+C/2„2 

' " kp 1/2 ^ /2 K 3/2H/2 (p) + f /2 -< /2 K 1/2+a2 (p)) . (4) 



3V2^r(l - C/2) V 2 

Here, p = pr d = p/p d and K u (p) is the modified Bessel function [|17|. A simi- 
lar Fourier transformation of the longitudinal structure function was performed by 
Sirovich, Smith, and Yakhot 0. []. Expanding eq. (§) for small p < 1 and ( > 
gives 

m = W- 1 ' 2 ^ r(3/2 + Q2) f 2-( \ 

L[P> 3s/2^,f r(l-(/2) P y + 2«l + P J' { ' 

i.e., we have a positive correction term to the expected asymptotic scaling E(p) oc 
p-C-i_ This correction signals the onset of an energy pileup around p d , see fig. 2. 
For large p>l the spectrum decays as E(p) oc p 1- ^ 2 exp (— p). Fig. 2 also shows a 
frequently used parametrization [|18|] for E(p), 



E(p) = cp <: exp(-p/p' d ), (6) 

where p' d is chosen in such a way that the r-space crossover r' d corresponding to eq. 
(^|) coincides with r d , for details see ref. |12|. This comparison emphasizes the energy 
pileup around p d described by the (modified) Sirovich-Smith- Yakhot formula eq. (|J), 
which can be considered to be an experimental spectrum summarizing the various 



simulations and experiments of ref. w and also those summarized in [14 



1 When the transcriptional error in eq. (20) of ref. H is corrected, the bottleneck pileup also 
shows up. 
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As already stated above, the energy pileup leads to a non monotonous local slope 
d In E(p) 



dlnp 



= -Up) - 1 



K(i + C)p 1/2 - c/2 - P 5/2 ~ c/2 ]K V2+(/2 (p) + (2 - C)p 3/2 ~ c/2 ^i/ 2+ c/2(p) m 

Cp 1/2 ^ /2 ^3/2 +C /2(p)+P 3/2 - f/2 ^l/2 + C/ 2 (P) 

For ( = 2/3 the maximum local slope is —1.448 (instead of —5/3) and occurs at 
Pmin ~ 0.85pd ~ (10t7) —1 . Fig. 3 shows ( p (p) together with ( r (r = 1/p) from eq. 
(0), demonstrating the strikingly different behavior of the local slopes in r- and in 
p-space. 

The energy pileup around pd has also been observed in further experiments Q 
(fitted by a correction term oc p 2 ^ 3 instead of our p 2 , cf. eq. @), in full numerical 
simulations [fj, |(| , and in a reduced wave vector set approximation (REWA) of the 
Navier-Stokes equations [f|. In ref. || a correction term +2(p/p pea k) 1 ' 8 was fitted to 
the data in nice agreement with our present result, +3p 2 /5 = +2.6(p/p pea k) 2 , where 
Ppeak is the point of maximum energy dissipation. Falkovich introduced the name 
"bottleneck phenomenon" for the energy pileup and predicts the correction term to 
be cx (p/pd) 4/3 /log(p d /p). 

We offer the following physical explanation (already given in of the bot- 
tleneck phenomenon: Consider the turbulent energy transfer downscale, T(p) ~ 
P u (p) I dpidP2u(pi)u(p 2 )S(p + Pi +P2), which does not depend on p in the inertial 
range due to Kolmogorov's structure equation fTLj . Assume that the amplitudes 
i^Pi), u{p 2 ) with pi, P2 > Pd > P are already damped by viscosity. Then the energy 
transfer T(p) would be reduced and stationarity could not be achieved unless u(p) 
increases. Because of the locality of the Navier-Stokes interaction in p-space, the 
effect is strongest around pd, leading to the energy pileup. Of course there is also 
viscous damping, but for p < pd < r/^ 1 it is smaller than the eddy viscosity T(p) 0. 
Borue and Orszag's simulations indeed show that the pileup starts in a region 
where T(p) is still constant. Above explanation rules out spectra of type (||). For 
an explanation of the bottleneck effect within the test field model we refer to ref. 
19|| , see also P0|| . For an analogous phenomenon in temperature spectra see J2T . 



Formally the bottleneck phenomenon reflects the relatively sharp crossover from 
r 2 -scaling (VSR) to r^-scaling (ISR) in the structure function (|1|). To illustrate this 
we transform the spectrum (^) back to r-space. This spectrum does not show the 
bottleneck phenomenon and the corresponding structure function 

D(r) = ^ Cr( ~ C ^ +1 (p' d r(C + 1) " (l + P' 2 d r 2 f +1) ' 2 sin ((C + 1) arctan^r))) 

(8) 

differs from the Batchelor parametrization ([!]) by its considerably smoother transi- 
tion (see the dotted curve in fig.l, showing a ratio of ~ 1.8 around r d -)- 

Our explanation suggests that the bottleneck effect potentially accompanies any 
sudden change in spectral strength, provided the wave vector amplitudes interact 
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nonlinearly and a conserved flux exists. We are consequently led to expect a similar 
effect at the infrared end of the scaling regime where the small-p modes are reduced 
in their spectral strength by the finite system size. 

Let us therefore consider the crossover between ISR and the large r saturation 
domain, where D{r) = 2(u 2 ) = 6w 2 rms becomes constant. Recall that L = 1/pl is 
the forcing scale. From experimental data fill , HI H3 we conclude that the second 
crossover at r = L is again well described by a Batchelor type transition 

D(r) = 2(u 2 )r 2 ■ (r 2 d + r 2 )" 1 ^ 2 ■ (L 2 + r 2 )"^ 2 , (9) 

see fig. 1. The general mechanism outlined above should equally well apply in this 
regime: The velocity amplitudes of the modes pi, P2 < Pl < P (or either of them) 
are reduced because of the finite size of the system. The mode u(p) again has to 
increase in order to guarantee a p-independent energy flux, now resulting in a steeper 
spectrum. 

Indeed, we find such a behavior for the spectrum corresponding to (^J). For 
r,j < r we derive the analytical result (for ( = 2/3) 



E(p) = <^(_E2M^ 



r(i/3) 



5_ 2 ( u 5 5 f £7 7^, 

g P i*2{ 6 , 2 , 2 , 4 J + m P i*2K 6 ' 2 ' 2' 4' 



71 

+ 2 



1 „ A 3 p\ 2 „ ,1 5 



(10) 



where p = p/pl and i-F 2 (a, b, c, z) denotes a generalized hypergeometric function 
17fl . The spectrum and the corresponding ( p (p) are shown in the left parts of figs. 



2 and 3, respectively. We find p max ~ 8pL and ((p ma x) ~ 0.77. Thus the deviations 
from classical scaling are again much larger than the discussed intermittency cor- 
rections. Note that our result agrees with theoretical [0, [TIJ and experimental hints 
(summarized in 0]) that the spectra are steeper for small p. 

We finally calculate the effective scaling exponent C^ e ^\Re\) that will be mea- 
sured in p-space simulations. Here we only consider the bottleneck phenomenon 
for large p, as in most numerical schemes the smallest wave vectors are forced and 
no p < pi are included. Let us express in terms of L and the Taylor- Reynolds 
number Re\ = \ui trms /u, where A = ui >rms / '{d\Ui) rms is the Taylor length. We have 
e = c t uf rms /L with c e = (6/6) 3//2 ~ 1, which is also known from grid turbulence 
experiments |22|. On the other hand, e = 15u(diUi) 2 ms pT| . Using these relations 
we finally get 77 = l^^c^LRe^' 2 or (for ( = 2/3) r d = (3&) 3/4 r? w 63LRe^ 3/2 . This 
connection between r d /L and Re\ allows us to calculate C^ e ^\Re\) as the average 

(^ f \Re x ) = 1 f ( p (p)dlnp, (11) 

\n(p min /p L ) J PL 



2 This crossover is probably nonuniversal. The important point here is simply the reduced 
spectral strength for small p, induced by the finite size, i.e., E(p) — > as p — » 0. 
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where p min is, as above, the wave vector of minimal ( p (p). The function (^ff\Re\) is 
shown in the inset of fig. 3. The deviations from the asymptotic value C,^ e ^\Re\ — ► 
oo) = ( are large. Assuming £ = 2/3, even for the largest Re\ = 200 achieved in 
numerical simulations |2| we have Cp 6 ^ ~ 0.58, which very well agrees with what is 
observed in numerical simulations 0. Impressive experimental confirmation of our 
prediction follows from recent measurements by Zocchi et al. ||. We include their 
data for (^ e ^\Re\) in our figure. 

Let us finally remark that our physical explanation of the bottleneck energy pile- 
ups is very general, it only assumes some inertial range with a constant energy flux 
in p-space. E.g., these conditions hold for surface or capillary waves p3| , where 
bottleneck phenomena are also expected ||, or for KS dynamics How bottle- 
neck phenomena manifest themselves in higher order moments and in power spectra 
remains a question for further research. 
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Figures 



Figure 1: Velocity structure function D(r), calculated from Benzi et al.'s data 
(dashed) [[[]], and its Batchelor fit ([]]) (solid) for ( = 2/3. Both curves are identical, 
which can even be seen in the enlargement of the crossover. The dotted line shows 
a structure function corresponding to the spectrum (||). The upper right part shows 
the saturation of the structure function (|9]), dashed. In the inset the longitudinal 
structure function Dl(t) is shown. The original data [I] are diamonds, the dashed 
line is a spline fit of these data, the solid line a fit of Batchelor type. Slight differences 
are seen. 



Figure 2: Experimental energy spectrum eq. (§) (solid) with ( = 2/3 and eq. (||) 
(dashed) without the energy pileup. In the left part the spectrum according to ( |T0| ) 
is shown. In the insets, the spectrum is enlarged around the energy pileups and 
compared to classical — 5/3-scaling. 



Figure 3: The local p-space scaling exponents ( p (p) from eq. (0) and eq. fllCf ) 
(solid), and the local r-space scaling exponent Cr( r Pd = 1/p) from eq. (0) and (|9p 
(dashed). The inset shows the averaged p-space scaling exponent C,^^\Re\) , see eq. 
(|TT|) (solid). Also shown is the local p-space scaling exponent ( p (pL(Re\)), dashed. 
We chose ( = 2/3 throughout. The dots are the experimental || (^ e ^\Re\). In || 
only (( e f^(Re) is given, so we calculated Re\ = cRe 1 ^ 2 with c = 0.17 chosen to give 
agreement for small Re\. 
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